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The existence of higher order entanglement in the stimulated and spontaneous Raman processes is established 
using the perturbative solutions of the Heisenberg equations of motion for various field modes that are obtained 
using the Sen-Mandal technique and a fully quantum mechanical Hamiltonian that describes the stimulated and 
spontaneous Raman processes. Specifically, the perturbative Sen-Mandal solutions are exploited here to show 
the signature of the higher order two-mode and multi-mode entanglement. In some special cases, we have also 
observed higher order entanglement in the partially spontaneous Raman processes. Further, it is shown that 
the depth of the nonclassicality indicators (parameters) can be manipulated by the specific choice of coupling 
constants, and it is observed that the depth of nonclassicality parameters increases with the order. 

PACS numbers: 03.67.Bg, 03.67.Mn, 42.50.-p 

Keywords: Higher order nonclassicality. Entanglement, Raman system. 


I. INTRODUCTION 

With the advent of quantum computation and communi¬ 
cation, entanglement has appeared as a very important re¬ 
source M. For example, its essential role in many pro¬ 
cesses, such as teleportation IH], dense coding ||2l, quantum 
information splitting IH], etc., are now well established. In 
short, entangled states are required to perform various impor¬ 
tant tasks related to quantum information processing. Entan¬ 
glement is produced in many physical systems and there ex¬ 
ists a large number of criteria for detection of entanglement 
(H and references therein). The first inseparability crite¬ 
rion was proposed by Peres |@] in 1996. Since then several 
inseparability inequalities have been reported for two mode 
and multi-mode states ilzHU- For the present study, we have 
mostly used higher order version of two criteria of Hillery and 
Zubairy Mill H^ . To be precise, we have used these crite¬ 
ria to investigate the existence of higher order entanglement 
in Raman processes, as depicted in Fig. [T] From Fig. [T] 
we can easily observe that the scheme illustrated here is es¬ 
sentially a sequential double Raman process S. Nonclas- 
sical properties of this system are studied since long (for a 
review see Ref. 1^ ). Initial studies on this system were re¬ 
stricted to the short-time approximation |2^ . However, 

recently nonclassical properties of this system have been in¬ 
vestigated by some of us ll^ l24l] using different approaches 
other than short-time approximation, but the possibility of ob¬ 
serving higher order entanglement is not investigated in any 
of the existing studies. Further, several applications of Ra¬ 
man processes have been reported in the recent past ll^ 1^ 
and higher order nonclassicality in different physical systems 
have also been reported experimentally 11271 - 1^ and theoreti¬ 
cally Keeping these facts in mind, present paper aims 

to investigate the possibility of higher order entanglement in 
the spontaneous, partially spontaneous and stimulated Raman 
processes and effect of the phase of the pump mode on the 
higher order entanglement. In what follows, Raman process is 
described as shown in Fig. [T]and a completely quantum me¬ 


chanical description of the system is used to obtain analytic 
expressions for the time evolution of the various filed modes 
involved in the process. The expressions are obtained using a 
perturbative method known as Sen-Mandal method 
Subsequently, the expressions obtained using this method and 
Hillery-Zubairy criteria are used to investigate the 

existence of multi-mode entanglement and higher order two¬ 
mode entanglement. Interestingly, the investigation has re¬ 
vealed the existence of multi-mode entanglement (which is 
essentially higher order as is witnessed via higher order cor¬ 
relation function) and higher order two-mode entanglement 
involving various modes present in Raman process. 



Figure 1: (Color online) Two-photon stimulated Raman scheme. The 
pump photon is converted into a Stokes photon and a phonon. The 
pump photon can also mix with a phonon to produce an anti-Stokes 
photon. 

Remaining part of the present paper is organized as follows. 
In Sectionini the Hamiltonian of stimulated Raman processes 
and its operator solution are briefly described. In Section HID 
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the solution is used to show the existence of higher order two¬ 
mode, three-mode, and four-mode entanglement and the effect 
of phase of the pump mode on the higher order entanglement. 
Finally, the paper is concluded in Section HVl 

II. MODEL HAMILTONIAN 

A completely quantum mechanical description of stimu¬ 
lated and spontaneous Raman processes described in Fig. [T] 
is given by the Hamiltonian 11^ 1^ 

H = uJaO^ a + uJbV^b + uJcC^ c +uJdd"^ d .... 

-I- g {ah^c^ + h.c.) + x {acd^ + h.c.) , 

where h.c. stands for the Hermitian conjugate. Throughout 
the present paper, we use h = 1. The annihilation (creation) 
operators a(a^), c (c^) , d{d^) correspond to the laser 

(pump) mode, Stokes mode, vibration (phonon) mode and 
anti-Stokes mode, respectively. They obey the well-known 
bosonic commutation relations. The frequencies uia, ^b, bJc 
and uJd correspond to the frequencies of pump mode a, Stokes 
mode b, vibration (phonon) mode c and anti-Stokes mode d, 
respectively. The parameters g and x are the Stokes and anti- 
Stokes coupling constants, respectively. Coupling constant g 
ix) denotes the strength of coupling between the Stokes (anti- 
Stokes) mode and the vibrational (phonon) mode and depends 
on the actual interaction mechanism. The dimension of g and 


X are that of frequency and consequently gt and xt are di¬ 
mensionless. Further, gt and xt are very small compared to 
unity. Further, we would like to note that in our present study, 
only one vibration (phonon) mode has been considered for 
the mathematical simplicity. In order to study the possibil¬ 
ity of the existence of higher order entanglement, we need 
simultaneous solutions of the following Heisenberg operator 
equations of motion for various field operators: 

a = —i {oJaO. + gbc -I- xcd) 
b = -i (ubb + gac^) 
c = —i (lOcC + gab^ -f 
d = -i {uJdd + xac) . 

The above set of equations (O is coupled nonlinear differen¬ 
tial equations of filed operator and are not exactly solvable in 
the closed analytical form under weak pump condition. How¬ 
ever, for the very strong pump, the operator a can be replaced 
by a c number and these equations (|2|i are exactly solvable in 
that case j^ . In order to solve these equations under weak 
pump approximation, we have used Sen-Mandal perturbative 
approach ll35l - [^ . The solutions obtained using this approach 
are more general than the one obtained for the same system 
using well-known short-time approximation. Details of the 
calculations are given in our previous papers ll35l - [38l] . Here 
we just note that under weak pump approximation, the solu¬ 
tions of Eq. (|2l) assume the following form: 


a{t) = /ia(0) -f /2&(0)c(0) -f f3C^{0)d{0) + /4a'l'(0)6(0)d(0) -f /5a(0)6(0)6'*'(0) 

+ /6a(0)ct(0)c(0) + /ra(0)ct(0)c(0) + /8a(0)dt(0)d(0), 
b{t) = gib{0) + p 2 a( 0 )c'l'( 0 ) -I- p 3 a^( 0 )(i'l'( 0 ) -I- g^c'^ (O)d(O) -f p56(0)c(0)c'l'(0) 

+ 56&(0)a(0)a'l'(0), ,,, 

c(f) = /iic(0)-f/i2a(0)6t(0)-f/i3at(0)d(0)-f/i4&^(0)ct(0)d(0)-f/i5c(0)a(0)at(0) ■ 

-f /i6c(0)6(0)&t(0) -b /i7c(0)dt(0)d(0) -b /i8c(0)at(0)a(0), 
d{t) = lid{0) + l2ai0)c{0) + l3a^i0)bH0) + Ub{0yi0) + Z5C^(0)c(0)d(0) 

-b (6a(0)at(0)d(0). 

The parameters /,, gi, hi and h are computed the initial boundary conditions. In order to obtain the solutions we use the 
boundary condition as at f = 0, in the first term of the Eq. Q. It is clear that /i(0) = 5i(0) = /ii(0) = Zi(0) = 1 and 
/i(0) = 5i(0) = hi{0) = li{0) = 0 (for z = 2, 3, 4, 5, 6, 7 and 8). Under these initial conditions the corresponding solutions 
for giit), hi{t) and li{t) are already reported in our earlier work Issl-fssl]. The same is included here as AppendixlAl 


III. HIGHER ORDER INTERMODAL ENTANGLEMENT 


In order to investigate the higher-order entanglement in 
spontaneous and stimulated Raman processes, we assume that 
all photon and phonon modes are initially coherent. In other 
words, the composite boson field consisting of photons and 
phonon is in an initial state which is product of coherent states. 
Therefore, the composite coherent state arises from the prod¬ 
uct of the coherent states Icti), |q; 2 ), jaa), and |a 4 ) which are 
the eigenkets of a, b, c and d respectively. Thus, the initial 


The solution Q, is valid up to the second orders in g and 
X- In what follows, we consider Awi = ojb + uJc — oJa and 
Aw 2 = oJa+uJc—iOd- The detunings Awi and Aa ;2 are usually 
very small. In the present work we have chosen |Aa;i| =0.1 
MHz and jAwal = 0.19 MHz. 
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composite state is 

IV'(O)) = |ai) (g) |a2) 0 lets) (g) |a4)- (4) 

It is clear that the initial state is separable. Now the field oper¬ 
ator a( 0 ) operating on such a composite coherent state gives 
rise to the complex eigenvalue ai. Hence we have, 

a(0)|V'(0)) = ailV'(O)), (5) 

where |ai p is the number of input photons in the pump mode. 
In a similar fashion, we can also describe three more com¬ 
plex amplitudes a 2 {t), a 3 (t) and 04 (f) corresponding to the 
Stokes, vibrational (phonon) and anti-Stokes field mode oper¬ 
ators b, c and d, respectively. It is clear that for a spontaneous 
process, the complex amplitudes except for the pump mode, 
are necessarily zero. Thus, in the spontaneous Raman process, 
= Q t = a 4 = 0 and ai ^ 0. For partially spontaneous 
process ll^ . the complex amplitude a\ and any one/two of 
the remaining three eigenvalues are nonzero while the other 
two/one complex amplitudes are/is zero. In the present inves¬ 
tigation, we consider that the eigenvalue corresponding to the 
pump mode is complex i.e., a\ = |ai| where is the 
phase angle, but the other eigenvalues (i.e., eigenvalues for the 
Stokes, vibrational (phonon) and anti-Stokes modes) are real. 

A. Higher order two mode entanglement 

In order to investigate the higher order two mode entangle¬ 
ment, we use two criteria due to Hillery and Zubairy lUlllO . 
The first criteria of Hillery and Zubairy is 

^ ^ q, (6) 


\ \ri2 / I |2(n-|-l) I |2(m-l) 

J = I/2I m\m\ai\ ' \a 2 \ 't 


and the second criterion is 


- 1 < 0 . ( 7 ) 


wherei and j are any two arbitrary operators and i,j € 
{a, b, c, d}\/i ^ j. Here m and n are the positive integers and 
the lowest possible values of m and nsae m = n = 1 , which 
corresponds to the normal (lowest) order intermodal entan¬ 
glement. A quantum state is said to be higher order entan¬ 
gled (bi-partite) if it is found to satisfy the equation (| 6 ]l and/or 
equation (|7]i for any choice of the integers m and n satisfying 
m + n > 3. From here onward we will refer to these criteria 
® and (|7]i as HZ-1 criterion and HZ-2 criterion, respectively. 
More specifically, a higher order entangled state is one which 
is witnessed via a higher order (order k > 2) correlation func¬ 
tion and as per this definition all multi-partite (multi-mode) 
entangled states are also higher order entangled. 

Before we proceed further, we note that these two criteria 
are only sufficient (not necessary) for detection of entangle¬ 
ment. Keeping this fact in mind, we have applied both of 
these two criteria to investigate the existence of higher order 
intermodal entanglement between various modes and have ob¬ 
served higher order intermodal entanglement in various situa¬ 
tions. In what follows, we have also investigated the possibil¬ 
ity of observing 3-mode and 4-mode entanglement. 

Let us first investigate the possibility of two mode entangle¬ 
ment in Raman process using HZ-1 and HZ-2 criteria. From 
Eqs. (O, (IHi, (|6]l and ([Til, we obtain the expression for the 
intermodal entanglement in pump and Stokes mode as 




In the similar manner, for the remaining cases, we obtain expressions for and j G {a, b, c, d}yi ^ j using HZ-1 

and HZ-2 criteria as follows 


rr^n.m 

rrifnm, 

c 


± 

± 

X 

X 

-f 


I52|^|a2|^^" (1 ± 2 to) |q;i|^ Ictal^-I-(1 ± 2n) |ai|^ |a 2 |^ 

w?!!? |ai|^ =F mn\a2'^ |a2|^” \ai'^ 

gig^mn a{a2a3 + gh^^mn 

(m — 1) (n — 1) -f (to — 1) \a 2 \^ -f (n — 1) |q! 3 |^| -f h 2 h\w?‘na\a 20 i\ 
1021^*-"“^^ + gigtrnn |q; 2|^*'"~^^ a2a|a4 ^2 |q; 3|^ -f to - 

h\^h 2 h 3 mn (to — 1) \ai\^ a 2 a^^a 4 + c.c. . 


(9) 
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rpn.m 

^a,c 

^a,c 


= |/2|^m|ai|^”|a3|^^™ (^m,\ai\^ ^ n\a^Y 

+ Ifsf \jn^ (1 ± 2 n) |ai|^ \a4\^ + (1 ± 2 m) jaal^ |a4f 

=F mn |ai|^ |q! 3 |^ + |a 4 |^| ± / 1 / 3 mn ai a 3 a 4 


/i 2 ^ 3 m^nQ;{^ ^2 |a 3 |^^™ Q ;4 + 72/3 |ai 012 |a 3 


ala\ 


+ fi^fl'm.naf' |a3|^^™”^^ a| |(n - 1) lasf + (m - 1) |ai|^ + ("»-iK"-i) | 

+ fif 2 hlh 3 mn {n - 1 ) q ;2 |Q!3|^"* 0:4 

+ fifshih^mn^m - 1) |ai|^"a 2 |a3|^^'"~^^ a§a 4 + c.c. 


( 10 ) 


77in,m 

^a,d 

T-^fn,m 

^a,d 


= |/3|^n|Q!i|^^” |a4|^'" |q; 4|^ T ™ |ai|^) 


( 11 ) 


J7in,m 

^c,d 

TTi/n.m 

^c,d 


= |/i 2 |" |ai|" |a 3 |"^”-'^ |a 4 |"'" + 1 ^ 21 " |a 4 |"'" |a 4 |" T mn |a 3 | 


( 12 ) 


'bd 

ifn,m 

b,d 


\g 2 \^ \ai\^ \a 2 \^^'^ |a 4 |^™ ± I\hmna\a 2 \a 2 \^^^ ^^ 041041 ^^™ + c.c. 


(13) 


Here we would like to note that once we obtain analytic ex¬ 
pressions for and in stimulated Raman process, it 

is straightforward to study the special cases: (i) spontaneous 
Raman process, where a 2 = 0:3 = 0:4 = 0, but ai ^ 0, 
and (ii) partially spontaneous Raman process, where cti 7 ^ 0 
and any one/two of the other three on {i = 2, 3, 4) is/are non¬ 
zero. It is clear from the Eqs. dSlflJll that for spontaneous Ra¬ 
man process Eqs. (ISlfTSTl reduces to zero. Hence for the spon¬ 
taneous Raman process, no signature of intermodal entangle¬ 
ment is observed. To investigate the possibility of higher order 
intermodal entanglement in the stimulated Raman process we 
have used y = g = 10^ Hz, |ai| = 10, |Qf 2 | = 8 , [asl = 0.01, 
1 0 : 4 1 = 1 . We have plotted the right hand side of dSl-fOl) in 

Eig. |2]and Eig. [3]for m = land n = 1, 2 and 3.We observed 
that HZ-1 criteria can detect the higher order intermodal en¬ 
tanglement in the stimulated Raman process for different val¬ 
ues of the phase angle or all phase angles of the input pump 
field (i.e., for </) = 0 , ^ and tt) for all the possible modes ex¬ 
cept pump-phonon (ac) and phonon-anti-Stokes (cd ) modes. 
It is interesting to note that higher order intermodal entangle¬ 


ment is observed in pump-Stokes mode, although in the low¬ 
est order it was not observed. Eurther, the figures show that 
the depths of the nonclassicality parameters and E'E’"^ 
increase with the order. Use of HZ-1 criteria also led to sim¬ 
ilar features in the partially spontaneous Raman process (not 
in figure). In other words, we observed signatures of inter¬ 
modal entanglement in all the cases except pump-phonon (ac) 
and phonon-anti-Stokes (cd) modes. As HZ-1 is only a suf¬ 
ficient (not necessary) criterion, it may have failed to witness 
entanglement, keeping this fact in mind, we have plotted the 
right hand side of Eq. ® fOT l using HZ-2 criteria (See Eig. 
O. It is interesting to note that HZ-2 criterion can detect the 
higher order intermodal entanglement in pump-phonon (ac) 
mode for phase angle (p = ^, which was not detected by 
HZ-1 criterion, in the stimulated and partial spontaneous Ra¬ 
man processes. However, we do not observe any signature of 
higher order intermodal entanglement for spontaneous Raman 
process. Thus, the stimulated Raman process provides a very 
nice example of a physical system which can produce higher 
order entanglement. 


J 


J 


B. Three mode entanglement 

There exists another alternative way to study the higher- 
order entanglement. To be precise, all multi-mode entan¬ 


glements are essentially higher-order entanglement. In other 
words, three mode entanglement always indicates higher or- 
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Figure 2: (Colour online) Fligher order intermodal entanglement in stimulated Raman process with |qi| = 10, |a2| = 8, |q3| = 0.01 
and 1 014 1 = 1 using HZ-1 criteron is shown for diffrent values of phase angle (0 = 0, and tt) in pump mode. Higher order intermodal 
entanglement is observed in (a) pump-Stokes mode for phase angle 0, (b) Stokes-vibration phonon mode for phase angle 0 = f ) (d) pump- 
anti-Stokes mode for phase angle 0, (f) Stokes-anti-Stokes mode for phase angle 0; and not observed in (c) pump-vibration phonon mode and 
(e) vibration phonon-anti-Stokes mode. In all the plots, the smooth line, dotted line and dash-dotted line are used for the m = 1 and n = 1, 2 
and 3, respectively. In (e), n = 2 and 3 are multiplied by 10® and 10®, respectively. While in all the remaining cases, n = 1 and 2 are shown 
1500 and 50 times, respectively. 


der entanglement. In order to investigate the three mode en- phonon mode respectively. Using equations (|3ll, (IHi and (fl4li 
tanglement, we use the following criterion we obtain 

EUc= {Na){m){N,) - \{abc)f < 0, (14) 

where (Na) , {N},) , and (Nc) are average value of the num¬ 
ber operators of the pump mode Stokes mode and vibration 

I 

E'a,b,a = (Na) im) (N,) - \{abc)f 

= \hf fs |ai|^ \a2f - \aif \az\^ - |q;i|'‘ - \ai\^ |a2|^) 

+ \h\^ (l«i|^ l«2|^ Ictal^ - 4 |q;2 |^ |Q!4|^ - 3 \a2\^ |a4|^ - 3 \aif \a2\‘^ Wi?') + 

- Iii/ija* |ai|^a2Q;3 + 2 / 1/3 cti \a2\^ ot^^aX + h2hl {2oL\aXaX 
-I- \ai'^ a\aXaX + ‘ 2 -a\aX\a 2 \^ a*^ + c.c. . 

For the spontaneous Raman process, Eq. (fTSl l reduces to 

K,.,c = -|/2|^l«l|^ 

which is clearly negative and thus indicate the existence of triparti te entanglement in the spontaneous Raman process 

To investigate the existence of three mode entanglement in the angle of the input pump field, i.e., for (j) = 0 (blue smooth 
stimulated Raman process, we plot the right hand side of the line), </ = § (red dotted line) and (j) = tt (green dash dotted 
equation (flSl l in Fig|4] for three different values of the phase line). The negative regions of the plots clearly illustrate the ex- 


(15) 


(16) 
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Figure 3: (Colour online) Higher order intermodal entanglement in stimulated Raman process is illustrated using HZ-2 criterion with different 
phase angle (0 = 0,^ and tt) in pump mode for |ai| = 10, |a 2 | = 8, |a 3 | = 0.01 and \ai\ = 1. Specifically, higher order intermodal 
entanglement is observed in (b) for Stokes-vibration phonon mode with phase angle (j) = f-, (c) for pump-vibration phonon mode with phase 
angle 0 = ^, and (f) for Stokes-anti-Stokes mode with phase angle ^. However, in (a), (d) and (e) higher order intermodal entanglement is 
not observed for pump-Stokes, pump-anti-Stokes and vibration phonon-anti-Stokes modes, respectively. The smooth, dotted and dash-dotted 
lines are used for m — 1 and n = 1, 2 and 3, respectively. Here, for (e) n = 2 and 3 are multiplied by 10® and 10®, respectively. For all the 
remaining cases, n = 1 and 2 are shown 1500 and 50 times, respectively. 


istence of tri-modal (higher order) entanglement. From Fig|4] 
we can clearly observe the signature of higher order entan¬ 
glement for different values of phase angle of the input pump 
field for stimulated and spontaneous Raman processes. 

C. Four mode entanglement 

In order to investigate the four mode entanglement we use 
the following criterion, which is similar to that of Li et al.’s 


three mode criterion lld^ : 

EUc,d = {m) (N,) (iVrf) - \{abcd)f < 0, (17) 

where a,b,c and d are arbitrary operators and the negative 
value of j ^ gives the signature of the higher order entan¬ 
glement. Now, we investigate the higher order entanglement 
i.e., the entanglement among the four modes of the stimulated 
Raman and spontaneous Raman processes and we obtain 


J 


Ea,b,c,d = \h\^ l«il^ (5 1^21^ \on\^ - |ai|^ - |ai|^ 

+ Ifsf \a2f \ai\^ (j \ai\^ - 4 \a4\^ - 3 \aif laif - 3 

- h\h2 \aif aia^a^ |a4|^ + 2f^f3a\ \a4f 04 + /z/g |a 2 |^ a 2 a§aj |q;4|^ 

+ ((4/3 - h 2 h^) |ai|^ ala 2 Ictal^ 0^4 + ^^2^3 {^alal^aX |a 4 |^ + |q;i|^ a{aXa\ \a 4 \^ ( 18 ) 

-F 2alat, |a 2 |^ 014 |a4|^ + 3a4a2 {asf aj |a!4|^) + fif2K^2alat,^aX'^ \a4\^ 

+ hXh2lilX aXaX^a4 + a\ \a2f \a4\^ q ; 4 ^ + 2f^f^lil^aX^ | q ; 2 |^ 

-F {h\h4 + ?>g\g 2 h\h 3 ) \aif \a4\^ a4 + c.c. 


In order to investigate the possibility of observing 4-mode 
entanglement in the Raman processes, in Fig. |5] we have 


plotted the variation of right hand side of Eq. (fTSl) with the 
rescaled time gt. Quite interestingly, for appropriate choice 
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Figure 4: (Colour online) The variation of three mode entanglement 
among pump, Stokes and vibration phonon modes (a) for stimu¬ 
lated Raman process using |ai| = 10, \a 2 \ = 8, las] = 0.01, 
|a 4 | = 1 (b) for spontaneous Raman process using |ai| = 10, 
\oL 2 \ = 1 03 1 = 1 04 1 = 0 with the smooth, dashed and dash-dotted 
lines corresponding to 0 = 0, and tt, respectively. 


Figure 5: (Colour online) Four mode entanglement among pump, 
Stokes, vibration phonon and anti-Stokes modes is depicted in (a) 
stimulated Raman process using |oi| = 10, |o2| = 8, |o3| = 0.01, 
1 04 1 = 1 (b) partial spontaneous Raman process using |oi| = 10, 
1 02 1 = 1 03 1 = 0. Flere, the smooth line, dashed line and dash-dotted 
line corresponds to ((> = 0, and tt, respectively. 


of the phase of the pump mode, 4 mode entanglement is ob¬ 
served in both stimulated Raman process and partially spon¬ 
taneous Raman process. 


IV. CONCLUSIONS 

Recently, nonclassical properties of the stimulated Raman 
process have been extensively studied by some of the present 
authors 1^ . In those studies intermodal entanglement 
in different modes of the stimulated Raman process was re¬ 
ported. Intermodal entanglement between Stokes mode and 
the vibration mode in the Raman processes was also reported 
by Kuznetsov ll43l] . However, higher order entanglement was 
not investigated. In the present paper higher order entangle¬ 
ment in stimulated Raman process is studied in detail and the 
observed higher order entanglement are illustrated through the 
negative regions of the FigsllJS] In Figs. |2] and [2 the ex¬ 
istence of higher order two-mode entanglement between vari¬ 
ous possible combinations of modes are illustrated using HZ-1 
criterion and HZ-2 criterion, respectively. Specifically, using 


HZ-1 criterion, we have observed the intermodal higher order 
entanglement for all the possible combinations of modes, ex¬ 
cept pump-phonon (ac) and phonon-anti-Stokes (cd) modes 
in stimulated Raman process (cf. Figl2]l and in partially spon¬ 
taneous Raman processes (not shown in figure). However, we 
found that HZ-2 criteria can detect the signature of higher 
order intermodal entanglement only in Stokes-phonon (6c), 
pump-phonon (ac) and Stokes-anti-Stokes {bd) modes in the 
stimulated and partially spontaneous Raman process Figl2 
but it is interesting to note that HZ-2 criteria can detect the 
higher order intermodal entanglement in pump-phonon (ac) 
mode whereas HZ-1 criteria fails to detect this. Thus, by 
combining the results, we have observed the existence of two¬ 
mode higher order entanglement in stimulated and partially 
spontaneous Raman possesses in all possible cases except in 
phonon-anti-Stokes (cd) modes. However, no signature of in¬ 
termodal entanglement is observed for the spontaneous Ra¬ 
man process. Another interesting point is that the present in¬ 
vestigation reveals the signature of higher order intermodal 
entanglement in pump-Stokes mode (a6) in stimulated Raman 
process, but intermodal entanglement in ab modes were not 
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observed in lowest order (cf. Fig. 2a, 3a and 4a of Ref. lH^l). 
As all the multi-partite (multi-mode) entanglement are essen¬ 
tially higher order entanglement, we investigated the possi¬ 
bility of observing three mode and four mode entanglements 
in Raman processes and found that tri-modal entanglement 
can be observed among pump, Stokes and vibration phonon 
mode (abc) in both stimulated and spontaneous Raman pro¬ 
cesses (cf. FiglDi, and it is also possible to observe entan¬ 
glement among four modes (pump, Stokes, vibration phonon 
and anti-Stokes) in stimulated and partially spontaneous Ra¬ 
man processes (see Fig|3. As recently many applications of 
multi-partite entanglement has been proposed, we hope that 
the present observation on the possibility of observing multi- 
mode entanglement in Raman process would be of help in re¬ 
alizing some of the recently proposed schemes that are based 
on multi-partite entanglement. Further, it is easy to experi¬ 
mentally realize Raman process and thus the results reported 
here can be experimentally verified using the available tech¬ 
nologies. 

Bosonic Hamiltonians similar to the one studied here fre¬ 
quently appear in quantum optical, opto-mechanical and 
atomic systems. Thus, the methodology adopted here may 
also be used in those systems to study the existence of non- 
classical states in normal and higher order entanglement in 
particular. Keeping this in mind, we conclude the present 
work with an expectation that this work would lead to a bunch 
of similar studies in other bosonic systems. 
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Appendix A: Parameters for the solutions in Eq. (3) 


gi = exp{-iu}bt), 
92 


93 = 


94 = 

95 

96 


ftn ioJbt 

r _ -i 1 

Acui ’ 

lpi{A(jJi—Auj2)t _ ll 

AcJ2(Aa;i—A ci;2) -I 

xge~“^i>* r iAiJit _ ll 
Aa;2Aa;i 


- 1 ]: 


_ AcJ; 


Xge \ i{Auji-\-Auj2)t _ i1 

)2(ALJi+Aa;2) l J 

r t -11 


xge [piAuit _ i] 

AUJ2A0JI L J ' 

Ji _ \piAu:it _ 1 ] _ ^9 te 

95- 


(A2) 


hi = exp(—zwcf) 
h2 = 


-iuict 

9^ 


Awi 


- 1 ] 


/iQ = — 


■x/c 'iujci 

_ |‘giAaj2i 


AiJo 


xge 


<^2 


/I 4 = 


Aa ;2 

XQe~ 


- 1 ] 

(Aaji +At>J2)* _ 


iAwi t 


hb 

he 

hr 


Aui 

g2^-iu:ct 

Aujf 


AcJi -\-AuJ2 

^i(A.uji+Auj2)t 


Acui 
iAuj2t 


Acui-|-A cJ2 


- 1 ] + 


Aa ;2 


ig te 


= —he 


^2^-iuict 

AujI 




hs = 


-\/2 o iujct 

^ ^ 3 rg*A^2t _ u _ 

Acoj L J 


Awi 

AuJ2 

AuJ2 


(A3) 


/i 

/2 

/3 

/4 

/s 

/6 

/? 

/s 


exp(-iwof), 

9e ^^-^Au:,t _ ^ 


Aloi 

All)' 

X.ge 


2 

— iuja t 


- 1 ] , 

-i(Aa;i-Aw2)t l 


AiiJi 


xge 


Au2 

g2g-i^a.t 

Aijf 


AijJi — Au)2 

■^(Aaji — AcAj2)t 


AcJs 

^-iAu;i t 


Acui — Au;2 


-iAoJit 


-i] + 


Aa;i 
2j.^ — iu)at 


ig te 


/s, 

^2^ — iuJat 
A ^ _ r iAuJ2t 

Awl L 
-/?■ 


Awi 

2j.„ — iuJat 


iX te 


Al02 


(Al) 


h 

h 


I 5 

h 


exp(-zwdf) 

o-iuidt 


Xe 


AuJ2 

xge 


^—iAuj2t 


- 1 ] 


7 _ AcJi (Acui —Aci;2 ) 

- —icojt r. 


^i{A(jJi — Auj2)t _ 


1 ] 


xge 




—iAuj2t 


- 1 ] 


u = 


Aa;2 Acui 

f -i(Auii+Aui2)t _ i ] 

AcJi (Acui+Aa;2) 


xge 


Auj2Auj 

ix^te~'‘‘^'^*' 


^ [e-*A-2t _ 1 ] 


Aw 2 


■Y/2p ioJdt 

+ ^ A 2 - ll 

Aw| L J 


= h 


(A4) 


1 ] 

















































9 


[1] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres and 
W. K. Wootters, Phys. Rev. Lett. 70, 1895 (1993). 

[2] C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 2881 
(1992). 

[3] C. Shukla and A. Pathak, Phys. Lett. A 377, 1337 (2013). 

[4] A. Pathak, Elements of Quantum Computation and Quantum 
Communication (CRC Press, Boca Raton, USA, 2013). 

[5] O. Gtihne, G. Toth, Phys. Rep. 474, 1 (2009). 

[6] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[7] L. M. Duan, G. Giedke, J. 1. Cirac and P. Zollar, Phys. Rev. Lett. 
84, 2722 (2000). 

[8] H. Hung and G. S. Agarwal, Phys. Rev. A, 49, 52 (1994). 

[9] J. Lee, M. S. Kim and H. Jeong, Phys. Rev. A, 62, 032305 

( 2000 ). 

[10] R. Simon, Phys. Rev. Lett. 84, 2726 (2000). 

[11] M. Hillery and M. S. Zubairy, Phys. Rev. Lett. 96, 050503 
(2006). 

[12] M. Hillery and M. S. Zubairy, Phys. Rev. A 74, 032333 (2006). 

[13] M. Hillery, H. T. Dung and H. Zheng, Phys. Rev. A 81, 062322 

( 2010 ). 

[14] G. S. Agarwal and A. Biswas, New J. Phys. 7, 211 (2005). 

[15] C. H. Raymond Ooi, Q. Sun, M. Suhail Zubairy and M. O. 
Scully, Phys. Rev. A 75, 013820 (2007). 

[16] C. H. Raymond Ooi, Phys. Rev. A 76, 013809 (2007); Eyob A. 
Sete, and C. H. Raymond Ooi, Phys. Rev. A 85, 063819 (2012). 

[17] A. Miranowicz et al., Phys. Rev. A 82, 013824 (2010). 

[18] T.-C. Wei and P. M. Goldbart, Phys. Rev. A 68, 042307 (2003) 

[19] A. Miranowicz, M. Piani, P. Horodecki, and R. Horodecki, 
Phys. Rev. A, 80, 052303 (2009). 

[20] A. Miranowicz and S. Kielich, Mdoern Nonlinear Optics Part 3, 
Ed. M. Evans and S. Kielich, Adv. Chem. Phys. 85, 531 (Wiley, 
New York, 1994). 

[21] P. Szlachetka, S. Kielich, J. Pefina and V. Perinova, J. Phys 
A: Math. Gen 12, 1921 (1979). 

[22] J. Perina, Quantum Statistics of Linear and Nonlinear Optical 
Phenomena (Kluwer, Dordrecht, 1991). 

[23] A. Pathak, J. Kfepelka and J. Pefina, Phys. Lett. A 377, 2692 
(2013). 

[24] B. Sen, S. K. Giri, S. Mandal, C. H. R. Ooi and A. Pathak, Phys. 


Rev. A 87, 022325 (2013). 

[25] P. Grangier, Nature 438, 749 (2005). 

[26] L. M. Duan, M. Lukin, J. 1. Cirac and P. Zoller, Nature 414, 413 
(2005). 

[27] A. Allevi, S. Olivares, and M. Bondani, Phys. Rev. A 85, 
063835 (2012). 

[28] A. Allevi, S. Olivares, and M. Bondani, Int. J. Quant. Info. 8, 
1241003 (2012). 

[29] M. Avenhaus, K. Laiho, M. V. Chekhova, and C. Silberhom, 
Phys. Rev. Lett 104, 063602 (2010). 

[30] E. Haas, J. Volz, R. Gehr, J. Reichel, and J. Esteve, Science 
344, 180 (2014). 

[31] S. K. Giri, B. Sen, C. H. R. Ooi and A. Pathak, Phys. Rev. A 
89, 033628 (2014). 

[32] K. Thapliyal, A. Pathak, B. Sen and J. Pefina, Phys. Rev. A 90, 
013808 (2014). 

[33] S. K. Giri, K. Thapliyal, B. Sen, and A. Pathak, 
arXiv: 1407.1780*^1 [quant-ph] 

[34] K. Thapliyal, A. Pathak, B. Sen and J. Pefina, Phys. Lett. A 
378, 3431 (2014). 

[35] B. Sen and S. Mandal, J. Mod. Opt. 52, 1789 (2005). 

[36] B. Sen, S. Mandal and J. Pefina, J. Phys. B:At. Mol. Opt. Phys. 
40, 1417 (2007). 

[37] B. Sen and S. Mandal, J. Mod. Opt. 55, 1697 (2008). 

[38] B. Sen, V. Pefinova, J. Pefina, A. Luks, and J. Kfepelka, J. 
Phys. B: At. Mol. Opt. Phys. 44, 105503 (2011). 

[39] P. Szlachetka, S. Kielich, J. Pefina and V. Pefinova, Opt. Acta 
27, 1609 (1980). 

[40] D. E. Walls, Z. Phys. 237, 224 (1970). 

[41] Same values of x, 9 and \ai\ are used in the entire paper (un¬ 
less otherwise specified). Eor spontaneous and partially spon¬ 
taneous Raman proceses these values of \ai \ are used for non¬ 
zero |ai|’s. 

[42] Z.-G. Li, S.-M. Pei, Z.-X. Wang and K. Wu, Phys. Rev. A 75, 
012311 (2007). 

[43] S. V. Kuznetsov, O. V. Man’ko and N. V. Tcherniega, J. Opt. B: 
Quantum Semiclass. Opt. 5, S503 (2003). 



